Search for: All records

Abstract Using very-long-baseline interferometry (VLBI) observations at (sub)millimeter wavelengths, the Event Horizon Telescope (EHT) currently achieves the finest angular resolution of any astronomical facility, necessary for imaging the horizon-scale structure around supermassive black holes. A significant calibration challenge for high-frequency VLBI stems from rapid variations in the atmospheric water vapor content above each telescope in the array, which induce corresponding fluctuations in the phase of the correlated signal that limit the coherent integration time and thus the achievable sensitivity. In this paper, we introduce a model that describes station-based phase corruptions jointly with a parameterization for the source structure. We adopt a Gaussian process (GP) prescription for the time evolution of these phase corruptions, which provides sufficient flexibility to capture even highly erratic phase behavior. The use of GPs permits the application of a Kalman filtering algorithm for numerical marginalization of these phase corruptions, which permits efficient exploration of the remaining parameter space. Our model also removes the need to specify an arbitrary “reference station” during calibration, instead establishing a global phase zero-point by enforcing the GPs at all stations to have fixed mean and finite variance. We validate our method using a real EHT observation of the blazar 3C 279, demonstrating that our approach yields calibration solutions that are consistent with those determined by the EHT Collaboration. The model presented here can be straightforwardly extended to incorporate frequency-dependent phase behavior, such as is relevant for the frequency phase transfer calibration technique. 

Abstract Very long baseline interferometry (VLBI) provides the highest-resolution images in astronomy. The sharpest resolution is nominally achieved at the highest frequencies, but as the observing frequency increases, so too does the atmospheric contribution to the system noise, degrading the sensitivity of the array and hampering detection. In this paper, we explore the limits of high-frequency VLBI observations usingngehtsim, a new tool for generating realistic synthetic data.ngehtsimuses detailed historical atmospheric models to simulate observing conditions, and it employs heuristic visibility detection criteria that emulate single- and multifrequency VLBI calibration strategies. We demonstrate the fidelity ofngehtsim’spredictions using a comparison with existing 230 GHz data taken by the Event Horizon Telescope (EHT), and we simulate the expected performance of EHT observations at 345 GHz. Though the EHT achieves a nearly 100% detection rate at 230 GHz, our simulations indicate that it should expect substantially poorer performance at 345 GHz; in particular, observations of M87* at 345 GHz are predicted to achieve detection rates of ≲20% that may preclude imaging. Increasing the array sensitivity through wider bandwidths and/or longer integration times—as enabled through, e.g., the simultaneous multifrequency upgrades envisioned for the next-generation EHT—can improve the 345 GHz prospects and yield detection levels that are comparable to those at 230 GHz. M87* and Sgr A* observations carried out in the atmospheric window around 460 GHz could expect to regularly achieve multiple detections on long baselines, but analogous observations at 690 and 875 GHz consistently obtain almost no detections at all. 

ABSTRACT We present MeqSilhouette v2.0 (MeqSv2), a fully polarimetric, time-and frequency-resolved synthetic data generation software for simulating millimetre (mm) wavelength very long baseline interferometry (VLBI) observations with heterogeneous arrays. Synthetic data are a critical component in understanding real observations, testing calibration and imaging algorithms, and predicting performance metrics of existing or proposed sites. MeqSv2 applies physics-based instrumental and atmospheric signal corruptions constrained by empirically derived site and station parameters to the data. The new version is capable of applying instrumental polarization effects and various other spectrally resolved effects using the Radio Interferometry Measurement Equation (RIME) formalism and produces synthetic data compatible with calibration pipelines designed to process real data. We demonstrate the various corruption capabilities of MeqSv2 using different arrays, with a focus on the effect of complex bandpass gains on closure quantities for the EHT at 230 GHz. We validate the frequency-dependent polarization leakage implementation by performing polarization self-calibration of synthetic EHT data using PolSolve. We also note the potential applications for cm-wavelength VLBI array analysis and design and future directions. 

Context.Sagittarius A* (Sgr A^*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. Aims.We investigated the variability and polarization properties of Sgr A^*using ALMA observations during the 2018 Event Horizon Telescope campaign. Methods.We analyzed high-cadence full-polarization light curves from ALMA at millimeter wavelengths, performed time-series analysis, and investigated the temporal behavior during an X-ray flare observed byChandraon 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. Results.We find low variability in total intensity (σ/μ < 10%), but significantly higher variability in linear and circular polarization (∼30% and ∼50%, respectively). A time-series analysis reveals red-noise variability, with power spectral densities between −2 and −3 across all Stokes parameters. Polarized intensity shows stable intra-day timescales, while total intensity exhibits more variable timescales, suggesting distinct emission regions, with polarization likely arising from a coherent structure. On April 24, a statistically significant inter-band delay in polarized intensity coincides with a near-simultaneous X-ray and millimeter peak that deviates from the typical delayed flare scenario. This event also features enhanced millimeter variability and coherent polarization loop evolution. The observed simultaneity challenges standard models of transient synchrotron emission with cooling delays, favoring instead a scenario of continuous energy injection in an optically thin region. Conclusions.Our results offer new constraints on the physical mechanisms driving variability in Sgr A^*, and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes. 

Abstract Event Horizon Telescope (EHT) images of the supermassive black hole M87* depict an asymmetric ring of emission. General relativistic magnetohydrodynamic (GRMHD) models of M87* and its accretion disk predict that the amplitude and location of the ring’s peak brightness asymmetry should fluctuate due to turbulence in the source plasma. We compare the observed distribution of brightness asymmetry amplitudes to the simulated distribution in GRMHD models, across varying black hole spina_*. We show that, for strongly magnetized (MAD) models, three epochs of EHT data marginally disfavor ∣a_*∣ ≲ 0.2. This is consistent with the Blandford–Znajek model for M87’s jet, which predicts that M87* should have nonzero spin. We show quantitatively how future observations could improve spin constraints and discuss how improved spin constraints could distinguish between differing jet-launching mechanisms and black hole growth scenarios. 

We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the significantly enhanced (u,v) coverage in the 2021 Event Horizon Telescope (EHT) observations. The integration of the 12−m Kitt Peak Telescope (USA) and NOEMA (France) stations into the array introduces two critical intermediate-length baselines to SMT (USA) and IRAM 30−m (Spain), providing sensitivity to emission structures at spatial scales of ∼250 μas and ∼2500 μas (∼ 0.02 pc and ∼ 0.02 pc). Without these new baselines, previous EHT observations of the source in 2017 and 2018 lacked the capability to constrain emission on large scales, where a “missing flux” of order ∼1 Jy is expected to reside. To probe these scales, we analyzed closure phases–robust against station-based gain calibration errors–and model the jet base emission using a simple Gaussian component offset from the compact ring emission at spatial separations > 100 μas. Our analysis revealed a Gaussian feature centered at (ΔRA ≈ 320 μas, ΔDec. ≈ 60 μ as), projected separation of ≈ 5500 AU, with an estimated flux density of only ∼60 mJy, implying that most of the missing flux identified in previous EHT studies had to originate from different, larger scales. Brighter emission at the relevant spatial scales is firmly ruled out, and the data do not favor more complex models. This component aligns with the inferred position of the large-scale jet and is therefore physically consistent with the emission of the jet base. While our findings point to detectable jet base emission at 230 GHz, the limited coverage provided by only two intermediate baselines limits our ability to robustly reconstruct its morphology. Consequently, we treated the recovered Gaussian as an upper limit on the jet base flux density. Future EHT observations with expanded intermediate baseline coverage will be essential to constrain the structure and nature of this component with higher precision. 

We present the first Event Horizon Telescope 1.3 mm observations of the supermassive binary black hole candidate OJ 287. The observations achieved an unprecedented angular resolution of 18 μas and reveal significant structural and polarization variability over just five days, marking the shortest timescale on which such changes have been directly imaged in this source. The inner jet exhibits a twisted ridgeline structure, with features displaying apparent superluminal motions up to about 22 c. The linear polarization maps reveal three main polarized features whose electric-vector position angles (EVPAs) change substantially over the time span of our observations, including a component with a radial polarization consistent with being produced by a recollimation shock. Most notably, we directly resolved two innermost jet components whose EVPAs rotate in opposite directions. The faster component, moving at 2.4 ± 0.9 μas/day (17.4 ± 6.5 c), exhibits counterclockwise EVPA swings of roughly 3.7° per day, while the slower component, with a proper motion of 1.4 ± 0.3 μas/day (10.2 ± 2.2 c), rotates clockwise at approximately 2.5° per day. Previous studies inferred helical magnetic fields in AGN jets from time-resolved or integrated polarization variability but lacked the angular resolution to directly image this effect. Our results provide spatially resolved evidence that a helical magnetic field threads the jet’s collimation and acceleration zone, ruling out models based on the superposition of unresolved components. Our analysis suggests that propagating shocks interact with a Kelvin–Helmholtz plasma instability, illuminating different phases of the helical magnetic field and producing the observed polarization spatial and temporal variability. Moreover, our model naturally accounts for the more rapid polarization rotation observed in the faster moving component. Our model predicts even more rapid swings in polarization, which could be tested with future observations featuring a more densely sampled time coverage. 

Aims.We investigated the polarization and Faraday properties of Messier 87 (M87) and seven other radio-loud active galactic nuclei (AGNs) atλ0.87 mm (345 GHz) using the Atacama Large Millimeter/submillimeter Array (ALMA). Our goal was to characterize the linear polarization (LP) fractions, measure Faraday rotation measures (RMs), and examine the magnetic field structures in the emission regions of these AGNs. Methods.We conducted full-polarization observations as part of the ALMA Band 7 very long baseline interferometry (VLBI) commissioning during the April 2021 Event Horizon Telescope (EHT) campaign. We analyzed the LP fractions and RMs to assess the nature of Faraday screens and magnetic fields in the submillimeter emission regions. Results.We find LP fractions between 1% and 17% and RMs exceeding 10⁵ rad m⁻², which are 1–2 orders of magnitude higher than typically observed at longer wavelengths (λ>3 mm). This suggests denser Faraday screens or stronger magnetic fields. Additionally, we present the first submillimeter polarized images of the M87 jet and the observed AGNs, revealing RM gradients and sign reversals in the M87 jet indicative of a kiloparsec-scale helical magnetic field structure. Conclusions.Our results provide essential constraints for calibrating, analyzing, and interpreting VLBI data from the EHT at 345 GHz, representing a critical step toward submillimeter VLBI imaging. 

We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87^*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT)'s 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measured the ring's ellipticity using five independent imaging methods, obtaining a consistent average value ofτ = 0.08_−0.02^+0.03with a position angle ofξ = 50.1_−7.6^+6.2 degrees. To interpret this measurement, we compared it to general relativistic magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including the thermal or nonthermal electron distribution function, spins, and ion-to-electron temperature ratios in both low- and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in nonthermal models and models with higher jet emission. These results indicate that the ellipticity measured from the M87^*emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.